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Abstract 

We consider a generalized chGUE based on a weak confining potential. We study the 
spectral correlations close to the origin in the thermodynamic limit. We show that for 
eigenvalues separated up to the mean level spacing the spectral correlations coincide with 
those of chGUE. Beyond this point, the spectrum is described by an oscillating number 
variance centered around a constant value. We argue that the origin of such a rigid 
spectrum is due to the breakdown of the translational invariance of the spectral kernel in 
the bulk of the spectrum. Finally, we compare our results with the ones obtained from 
a critical chGUE recently reported in the literature. We conclude that our generalized 
chGUE does not belong to the same class of universality as the above mentioned model. 
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1 Introduction 



Critical statistics has been an intense subject of study in recent years 0, |2^, 

T6| , ^. So far, two kinds of models have been proposed to describe those critical 



correlations. In one of them, deviations from Wigner-Dyson statistics are obtained by 
adding a symmetry breaking term to the Gaussian Unitary Ensemble (GUE) pO| , T0[ . 
The model is solved by mapping it to a non interacting Fermi gas of eigenvalues. The 
second one |^ makes use of soft confining potentials. It is solved exactly by means of 
q-orthogonal polynomials. 

Universality in critical statistics have been conjectured |^| due to the fact both models 
share the same kernel in the thermodynamic limit and when the deviations from GUE are 
small. However, the origin of the critical kernel is different in both cases. In models based 
in soft confining potentials the critical kernel is obtained from a nontrivial unfolding caused 
by the strong fluctuations of the spectral density. In models with a explicit breaking 
symmetry term, deviations from Wigner-Dyson statistics arise because the long range 
interactions among eigenvalues are suppressed |TT|. Although progress have been recently 
reported , the universality class associated with critical statistics can be still considered 
an unresolved problem. 

Recently a critical chiral Gaussian Unitary Ensemble (chGUE) of the first 
class (addition of a symmetry breaking term to the chGUE) was proposed in order to 
describe the spectral correlations of the QCD Dirac operator beyond the Thouless energy 
1^. It was found that in the bulk of the spectrum and for small deviations from the 
Wigner-Dyson statistics, the spectral kernel coincided with the one conjectured ||2T]] to be 
universal for critical statistics. In this letter we shall study the effect of the hard edge 
(the ensemble is defined on the positive real axis only) on the spectral correlations of a 
chGUE with a weak confining potential. 

We will proceed as follows. First, we propose a random matrix ensemble defined on the 
positive real line with a non-polynomial potential which is soft confining in the bulk of the 
spectrum and Gaussian close to the origin. Then, we compute the spectral kernel in the 
semiclassical approximation. Finally, we compare our model with the above mentioned 
critical chGUE 

Properties of chGUE with weakly confining potentials have already been discussed in 
the physics literature ^, |T3|, but attention was focused in the bulk of the spec- 
trum. The effect of the hard edge in the critical spectral correlations and its impact on 
universality remains an open question. 
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2 Definition of the Model 



In this section we introduce the model to be studied and argue the need to unfold the 
spectrum. Finally, we compute the mean spectral density needed for such unfolding by 
using the Dyson's mean field equation. 

We consider a N x N complex hermitian matrix ensemble H with block structure, 

and probabihty distribution given by, 

P(C) oc e-^(^^') (2) 

where C is a N/2 x N/2 hermitian matrix. In terms of the eigenvalues of H, the joint 
distribution is given by, 

N 

P{xi . . . xat) OC n a;.e-^("') [] 1^' " ^jT (3) 

1=1 l<i<j<N 



V{xi) = — arcsinh^(a;j) (4) 
7 

where Xj are the eigenvalues of H. 

Since V{xi) in (|^) is proportional to xf for -C 1 we expect to recover the chGUE 
kernel ||25[ in this limit. For Xi ^ 1, the potential V{xi) oc log^(xj) fail to keep the 
eigenvalues confined and deviations from the chGUE may be relevant . 

If the considered interval were the whole real line, the orthogonal polynomials associ- 
ated with (^) would be the 1/q- Hermite polynomials hn{x; q) 0, ^ with 7 = log(l/g). 
Unfortunately, for the positive real axes we do not know any set of polynomials orthogonal 
with respect to the measure (0) with the potential (^). Thus, in order to compute the 
mean spectral density necessary to unfold the spectrum we shall use the Dyson's mean 
field equation. 

The joint distribution @) can be written as a statistical distribution of a one di- 
mensional system of N particles at temperature "1/T = 2" with a pairwise logarithmic 
interaction and a one particle potential given by that maintains the system confined. 

P(xi,...x^) oce-^("^'-'"'^)/^ (5) 

where 

N N 

F{xi, ...,xn) = Y. - 2 ^ log -x^^\- log(x,) (6) 

i=l i^j 
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We want to perform a mean field theory analysis of the above one dimensional system. We 
assume that in the large limit the above system has a continuous macroscopic density 
given by, 

JV 

p{e) = Y.^{e-e.) (7) 

i=l 

Plugging p(e) into (Q) and assuming that the density is non zero only in the interval 
< X < D we can express F as a functional of the spectral density. 

FW)] = r p{e)V{e)de - T T p{e)p{e') log \e - e'\dede' (8) 
JO Jo Jo 

The mean spectral density pmf is defined as the density that minimizes the above 
functional, namely, 6F/Sp = 0, that implies 

r de'p{e',,p)log\e-e'\ = Vie)+c (9) 
Jo 

where c is a constant due to the normalization constraint. The general solution of the 
above equation (usually called Dyson equation) is given by [1T2| , 



1 ^ fD dV/dt tdt , , 

where e+ = e + iO. D is found from the normalization condition, 

PMF{e)de = N (11) 

Jo 

Now, the task is to compute pmf for the potential V{e) = iarcsinh^(e), 

f\ 2 rjr^ 20 /"^ arcsinh(t) 1 tdt 

77r^ Jo Vl + i V -D^ -t^t^ - €% 

This integral can be performed by changing the contour of integration in a sum of two 
pieces, A = Ai + A2 where Ai is the the negative imaginary axis and A2 is the interval 
[D,oo]. Since we are interested only in the real part of (|T2p, A2 does not contribute to 
the integral . Thus, (0) can be written as. 



p^,^(e) = —VD^ - e2- / = (13) 

The above integral can be performed by means of a change of variables, the final result 
being. 



1 1 VL)2 _ g2 

Pmf (e) = — r^— arctan ( 14) 

71^7 V e + 1 V e + 1 
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From the normalization condition we find that D = sinh(7r7A^), therefore, the mean 
spectral density for ^ oo is given by, 

Pmf e = - — T== 15 
27 V e + 1 

As expected, Pmf{^) has the right limiting values, it is a constant for e <^ 1 (as in the 
chGUE case) and for e ^ 1 is proportional to 1/e, as for the random matrix ensemble 



with soft confining potentials discussed in |T6|, |2^, ^. The above spectral density will 
be used in the next section to unfold the spectrum. 

This unfolding allows us to work in units in which the mean level spacing is equal 
to one. We recall that, in this context, random matrix theories only reproduce spectral 
correlations around the average spectral density. 

We remark that the above mean spectral density is an approximate formula capable 
of giving only the smooth part of the spectral density. The exact mean spectral density 
has oscillations which are out of reach of the mean field formalism above used. Therefore, 
the mean spectral density (O) is only valid approximation if these fluctuation are small 



enough |TB[. In our model this situation corresponds with 7 <^ 1. For 7^1 the exact 
spectral density is a rapidly oscillating function. Hence, it is not possible to define a 
meaningful mean spectral density out of it |]16|, 



As the mean spectral density is not constant, the rescaling procedure is not trivial 
TT[| . The variable x in terms of which the spectral density becomes a constant, is the 



integrated mean spectral density. 



rE 

X = I pMF{(^)de (16) 







where Pmf{^) is the mean spectral density previously found. We shall see in the next 
section that this nontrivial unfolding is the main ingredient to get a non translational 
invariant kernel in the bulk of the spectrum. 

3 Calculation of the spectral kernel 

In this section we compute the spectral kernel in the semiclassical approximation. The 
semiclassical approximation in the GUE consists in substituting the wave function (or- 
thogonal polynomials times e"^'-^-*) appearing in the spectral kernel, after the Christoffel- 
Darboux formula is applied, by their WKB approximation. Because of the hard edge at 
a; = we cannot simply do a WKB approximation by replacing the wave functions by 
plane waves, but instead have to use Bessel functions. The kernel associated with (^) can 
be written in terms of the wave function as follows, 

^/j2N+liE)^lj2NiE')-^/j2N+liE')^2NiE) , V^2iV+l(^0^27v(^) + ^2iV+l(^0^2jv(^) 

^ ^ = n{E - E') + n{E + E') ^^^^ 
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Figure 1: We compare the analytical value of T2(x, 6.53) = \K{x, 6.53) p with the numeri- 
cal one. K{x, y) is given by (^) with 7 = 0.25. The numerical integration was performed 
by using the Metropolis algorithm for = 100 "particles" in the confining potential @. 
We have used the first 9-10'' sweeps to "warm up" the system and taken the average over 
the next 9-10^. We repeated the process four times. The final result is the the average 
of the four trials. We explicitly checked the agreement between numerical and analytical 
results for y < 20. 



where ip2N{,E) and ip2N+i{,E) are the even and odd large limit of the wave functions 
associated with the potential (^. In the semiclassical approximation those functions are 
given by 0, 



1p2N+l{E) 



s{E)Jo{ns{E)) 



;i8) 



s{E)MnsiE)) 



where Jq and Ji are Bessel functions and s{E) is defined by, pmf = ds/dE, with pmf 
the mean spectral density computed in the last section. It is clear ||^ that the above 
semiclassical expressions for the wave functions are correct for polynomial increasing po- 
tential. For soft confining potentials, according to 0, ^ these expressions are valid if 
the mean-field spectral density used to unfold the spectrum is close to the exact mean 
spectral density. In our model this happens whether 7 <C 1. Othe argument supporting 
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(|l^) comes from the asymptotic form of the orthogonal polynomials associated with a po- 
tential asymptotically proportional to log^(-E'). This problem have already been discussed 
in the literature 0. They found that for E 1, '^2n{E) oc cos(log(-E')/7). This result 
coincides with (^) in the limit considered. For E <^ 1, s{E) oc I/7 and we recover the 
chGUE result. As an additional check, we evaluate T2{x,y) = by numerical 

integration of the the joint distribution (^. Figure 1 shows that the agreement between 
numerical and analytical results is excellent. We recall that for polynomial- like increasing 
potentials the mean spectral density is a constant proportional to S> 1, s{E) is linear 
for E ^ and we recover the kernel of the chGUE |^ . 



Once we know the form of the kernel we can unfold the spectrum by using (|T6D and 

(0) Hill. 



E 2 1 

de = X (19) 

-fVTT^ 



E = sinh(a;7/2) (20) 

where for convenience we have set 7 = 7/4. The kernel in terms of the new, unfolded 
variables is given by, 

8 ^ smh((a; + ?/)7/4) cosh((x — y)7/4) 

Ji{Tcx)Jo{ny) - Ji{TTy)Jo{7Tx) ^ ^^^^ 



sinh((x — y)^/4:) cosh((a; + y)^ /A) 
As expected, for 7 — we recover the chGUE kernel. This kernel is already in a suitable 



form for comparison with the one previously found in |[T0[| . 



i^[x,y) ^^xy[ sinh((x + y)7/4) sinh((a; - y)7/4) ^ 

Even though both kernels reproduce the chGUE kernel for 7 — they are essentially 
different in the bulk of the spectrum, (p^) is translational invariant in the bulk of the 



spectrum, unlike (0), which is not. The origin of such non traslationally invariant kernel 
is due to the nontrivial unfolding induced by the mean spectral density. This unfolding 
prevents from vanishing the second term of the right hand of (|^) in the bulk of the 
spectrum. 

In the next section, we shall study the effect of the non-translational invariance of 
the kernel in the spectral correlations involving many levels by computing the number 
variance. 
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Figure 2: Number variance for 7 = log(l/g) = 0.1. It is linear for L ^ 1 in the critical 



cliGUE of [10|. In our generalized chGUE is almost constant in the bulk of the spectrum. 



4 Discussion of results 



In this section we shall see, by computing the number variance, that the spectrum of 
our model is more rigid than the chGUE one and essentially different from the models 
describing critical statistics. 

In order to observe deviations from chGUE prediction we are going to study long-range 
correlations of eigenvalues by studying the number variance in an interval [0, s\. 
The number variance is a statistical quantity which gives a quantitative description of the 
stiffness of the spectrum. The number variance is obtained by integrating the two-point 
correlation function including the self-correlations 

T?{L) = Tdx ("^ dy[6{x - y){p{x)) + R^ix.y)]. (23) 
Jo Jo 

For the Wigner-Dyson statistics the number variance is proportional to \og{L). Such 
weakly increasing number variance is not surprising as the eigenvalues repulsion produces 
a highly rigid spectrum. For the Poisson statistics the number variance is equal to L as 
expected for eigenvalues which are not correlated. Finally, a number variance proportional 
to xL (for L ^ 1 and x ^ 1) is a signature of critical statistics 0, [l^, |2^, |T^. The 
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Figure 3: 7 = log(l/g) = 0.3. A constant number variance is observed in the generalized 
chGUE when L ^ 1. That means the spectrum is even more rigid than the chGUE one. 



slope X is directly related to the multifractality of the wave functions of a disordered 
system at a delocalization-localization transition. 

A linear number variance for L ^ 1 with a slope x 1 ^^^s found in the generalized 



chGUE ||10[- However, as it can be observed in Figure 2 and 3, the number variance of 
our model is almost constant for L ^ 1. The oscillating behavior around a constant value 
is partially due to the self-interactions coming from the first term of the number variance. 

Apparently, this result is surprising because random matrix ensembles with broken 
time invariance based on potentials behaving as log^(x) asymptotically are supposed to 
have a linear number variance with slope x iii the bulk of the spectrum, which is a 
signature of critical statistics |2^. In principle, one may think that the presence of a 
hard edge at x = in our model does not affect the spectral properties in the bulk of the 
spectrum. We argue that this is not the case. 

The hard edge, combined with the soft-confining nature of the potential breaks up the 
translational invariance of the kernel (|2T| ) even in the bulk of the spectrum. In the bulk , 
the cluster function associated with the kernel (|21|)is given by. 



Y2{x,y) oc [ 



+ 



sin2(7r(a; + y)) _^ ^_ sm{iT{x + y)) 



sin(7r(x — y)) 



sin^(7r(x — y)) 

sinh^((x — y)7/4) ' cosli^ {{x — y)^ / A) ' cosh((a; — ?/)7/4) sinh((a; — y)7/4) 



](24) 



9 



where we have used the asymptotic expression of the Bessel functions. By performing 
elementary integrations, we observe that the leading contribution to the number variance 
for L ^ 1 coming from the first term of Y2 (translational invariant part) is xL where x 
is a function of 7 only. On the other hand, the leading contribution of the second term 
of Y2 (non-translational invariant part) to the number variance is ~xL- Therefore, both 
contributions cancel each other and we are left with a oscillating (around a constant value 
depending on 7) number variance coming from the third term of the cluster function. We 
point out that the above cancelation is mainly due to the non-trivial unfolding used. 

Roughly speaking, weak increasing potentials fail to keep the eigenvalues confined. As 
a consequence, the mean spectral density is in general a strongly oscillating function even 
in the thermodynamic limit. If 7 -C 1 the deviations from chGUE are small and we can 
still define a relevant smooth mean spectral density by using the mean field formalism 
0. The unfolding procedure using this mean spectral density breaks the translational 
invariance of the kernel in the bulk of the spectrum. This breaking of the translational 
symmetry produces a spectrum highly correlated and essentially different from the one 
reported in ||10| . 

We would like to mention that a similar result to the one found in this letter has 
been reported by Canali and Kravtsov [|13|, |lT| . They studied the spectral properties 
of a generalized GUE based on a weak confining potential with a log^(x) asymptotic as 
well. They noticed that in the bulk of the spectrum for ^ 00 and 7 << 1, the cluster 
function Y2{x,y) of that ensemble has strong correlations not only when x ~ y, but also 
when X ~ —y. The total cluster function is given by the following non translational 
invariant relation, 



Y2{x y) = 7^ sin^(7r(x - y)) cosh(x7/2) cosh(y7/2) ^^5) 
167r2 cosh^((x + ?/)7/4) sinh^((x — ?/)7/4) 



They showed that, due to this 'ghost' peak, the number variance depends on the 
interval in which it is calculated. 

If the interval is not symmetric with respect to the origin ([0, s] for instance), the 
system does not feel the strong (non-translational invariant) correlations at x = ~y. 
Then, the number variance goes asymptotically like x^ and the model is supposed to 
describe critical correlations. However, if the interval is symmetric with respect to the 
origin , the peak at a; = —y of the two point function Y2{x, y) has to be taken into account 
as well. This contribution drives the asymptotic form of the number variance to a constant 
value [^, [1^, in agreement with the results obtained in this letter. 



It is straightforward to connect the number variance in the bulk of the spectrum of 
27| , |13[ with the one studied in this letter. The asymptotic form of the number variance 



in a interval [0, s] associated with the first two terms of (^41) corresponds to the number 
variance in the interval [— s/2, s/2] of the above mentioned critical GUE. By changing 
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variables u = —x,v = y in the expression for the number variance of our model we recover 
the expression obtained by Canali and Kravtsov. The third term of produces the 
oscillating behavior observed only in our generalized chGUE. 

To sum up, due to the non translational invariance of the kernel contributions coming 
from the points x ~ —y have to be taken into account. These contributions make the 
linear term in the number variance vanish. 



5 Conclusions 

In this letter we have studied the effect of a hard edge in the spectral correlations of a 
chiral random matrix ensemble with a soft confining potential. We showed that beyond the 
Thouless energy the spectrum is characterized by an oscillating number variance around 
a constant value. The spectrum is even more correlated than the chGUE one. 

The linear term of the number variance characterizing critical statistics vanishes due 
to the non-translational invariance of the spectral kernel in the bulk of the spectrum. 
Thus, the generalized chGUE studied in this letter and the ensemble of (in which a 
linear number variance was found to be proportional to xL for L >> 1) belong to different 
universality classes [|I|]. 
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